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Exact diagonalizations of two dimensional small t-J clusters reveal dominant 
hole-hole correlations at distance y/2 (ie between holes on next nearest neighbor 
sites). A new form of singlet d x 2_ y 2 pair operator is proposed to account for the 
spatial extention of the two hole-bound pair beyond nearest neighbor sites. 
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Some time ago it has been proposed that the exchange of spin fluctuations could lead to 
superconductivity in the singlet d x i_ y i channel [I]]. Such mechanisms were recently discussed 
in connection with the new copper oxyde superconductors. Occurence of d-wave pairing in 
strongly correlated models have been further supported by quantum Monte Carlo calcula- 
tions |2|, weak coupling approches || and recent finite size scaling analysis ||. The issue of 
pairing between holes moving in an antiferromagnetic fluctuating background is of central 
importance in the field of high-T c superconductivity. Theoretical predictions f| for NMR in 
a d-wave superconductor seem to be consistent with the measured data ||. Also consistent 
with d-wave pairing let us briefly mention recent microwave experiments |7| showing a linear 
temperature behavior of the penetration depth and recent photoemission data || revealing 
nodes in the gap function. 

In standard notations the t-J Hamiltonian on a two dimensional (2D) reads 

H = J Y1{S\ • S i+ g- -THm+g) - t J2(4,a c i+?,<r + h.c), (0.1) 

where Si and c- lj(7 stand for the localized spins and the itinerant holes respectively. 

Numerical calculations of the binding energy between two holes on 4x4 @-12], V4~8x a/18 



HJ, V20 x V20 [|Tj] and y26 x y26 [|J square clusters suggested that binding occurs for 
realistic values of the antiferromagnetic exchange interaction (0.3 < J/t < 0.6) fL5| . The 
bound state is found to be a spin singlet with d x 2_ y 2 internal symmetry. As revealed by a 
quasiparticle pole in the dynamical pair correlation function |T7],rj] the bound pair exhibits 
a coherent motion characteristic of a unique (bosonic) particle of charge +2e which leads to 
flux quantization with a periodicity in flux of $o/2 

The hole-hole correlations C(R) in the two hole-ground state (GS) of the \/26 x \/26 
cluster are shown in Fig. 1 in good agreement with previous data on the 4x4 cluster 

|I(],|l2|] and with approximate results on the same cluster obtained from a restricted Hilbert 
space analysis |Tj|. Over a broad region, 0.2 < J/t < 0.9, the correlation at distance y/2 



dominates. Strong correlations are also found at distance 1 (i.e. for holes sitting on nearest 
neighbor sites) which become in turn dominant for large J/t (unphysical). If the two holes 
do form a bound state, as it is suggested by a finite size scaling analysis 0, we expect the 



correlation C(R) at small distance R as well as the mean- hole separation J3J to remain finite 
in the thermodynamic limit. 

Our purpose here is to define a proper pair operator which properly accounts for the 
internal structure of the pair, with strong density-density correlations at distance 1 and \/2. 
A generalized form of a pair operator can be defined as 
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where and ^c^ -2 ar e the Neel and two hole GS respectively on the N site cluster and 
Ar creates a pair of holes on two sites at a distance R. The sum can be restricted to a small 
set of lattice distances determined by the actual size of the pair. Fig. 1 suggests that, in 
practice, it is sufficient to include R — 1 and R = y2. The prefactors in (p.2|) fl§| are chosen 
in order to maximize the quasiparticle (QP) weight of the pair pO|fl, 
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From Eq. ( p.3|) it is clear that the QP weight increases with the number of pair operator A# 
included in the sum. Hence adding a contribution at R — v2 could significantly increase 
the weights Z^h calculated in Ref. |f] with the conventional d x i_ y i BCS spin singlet operator 
for nearest neighbor sites Ai = J2i c i,T( c i+x,j + c i-x,j — c i+y,| — c i-y,l)- 

Unlike for R=l, the pair operator at separation x + y can not be written in a simple 
BCS form because of the d x 2_ y 2 character of the GS wavefunction. Indeed, simple symmetry 
arguments imply that the pair operator must be odd under reflections with respect to the 
x = ±y directions. Let us define A^ as 

A^ = XX^i+y ~~ ^i+x) ' ^M+x+y — (S'i-x — ^i+y) ' ^i.i-x+y}) (0-4) 
i 

where = \c\^{a y a) a(T ic^ a i is the regular spin triplet pair operator This expression is 



very similar to the pair function proposed by Balatsky and Bonca |2^] for the ID t-J model 
in the context of odd-gap pairing. However, we stress that (p.4j) is, as far as symmetry 
properties are concerned, entirely similar to A x (spin singlet, d x 2_ y 2 spatial symmetry) as 
schematically depicted in Fig. 2 whereas the odd-gap singlet operator of |22j is odd under 
space inversion, r — > — r (p-wave character). 



In summary, the hole-hole density correlation function calculated exactly on a ^/2Q x \/26 
cluster shows further evidence in favor of a two hole d x 2_ y 2 bound state in the fluctuating 
antiferromagnetic background. In order to include the strong tendancy of the holes to stay 
together across the diagonal of a plaquette a new form of pair operator is introduced which 
couples the hole pair to the neighboring spins. Although its form is very similar to pair field 
operators introduced in the context of odd-gap pairing its symmetry properties are identical 
to the symmetry properties of the conventional BCS pair operator at distance 1. 
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FIGURE CAPTIONS 

Figure 1 

Hole density correlations for 2 holes on a \^26 x \^26 cluster. The symbols corresponding to 
the various hole- hole separations are shown on the plot. 
Figure 2 

Schematic picture of the d-wave pair operator around a lattice site i. The arrows 
correspond to the triplet pair operators (oriented) and the signs ± in parenthesis correspond 
to the phase factors of the spin operators in ( p.4|) . 
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